Computational study of microwave oscillations in absence of external field in 
nonstandard spin valves in the diffusive transport limit 



o 
o 

(N 



o 

(N 



o 
o 



> 

oo 
m 
cn 

o 

On 
O 



X 



E. Jaromirska, 1 P. Balaz, 2 L. Lopez Diaz, 1 and J. Barnas 2 ' 3 

1 Departamento de Fisica Aplicada, Universidad de Salamanca, Spain 
2 Department of Physics, Adam Mickiewicz University, Umultowska 85, 61-614 Poznan, Poland 
1 Institute of Molecular Physics, Polish Academy of Sciences Smoluchowskiego 17, 60-179 Poznan, Poland 

(Dated: July 21, 2009) 

An anomalous (inverse) spin accumulation in the nonmagnetic spacer may build up when the spin 
valve consists of magnetic films having different spin asymmetries. This leads to wavy-like depen- 
dence of spin-transfer torque on the angle between magnetizations, as predicted by spin-dependent 
diffusive transport model, and also confirmed experimentally. Making use of these predictions, we 
have numerically studied the magnetization dynamics in presence of such a wavy torque in Co(8 
nm)/Cu(10 nm)/Py(8 nm) nanopillar, considering geometry with extended and etched Co layer. 
In both cases we specify conditions for the out-of-plane precession to appear in absence of external 
magnetic field and neglecting thermal fluctuations. We prove the assumption of wavy-like torque 
angular dependence to be fully consistent with experimental observations. We also show that some 
features reported experimentally, like nonlinear slope of frequency vs. current behavior, are beyond 
the applicability range of macrospin approximation and can be explained only by means of full 
micromagnetic analysis. 

PACS numbers: 67.30.hj,75.60.Jk,75.70.Cn,78.20.Bh 



I. INTRODUCTION 

The concept of spin transfer was introduced in pio- 
neering works by Slonczewskii and Bcrgcr- 2 . They have 
shown that spin polarized current can exert torque on 
a thin magnetic film due to transfer of spin angular mo- 
mentum, influencing the magnetic state of the layer. The 
spin-transfer torque (STT) originates from spin asym- 
metries of the two independent transport channels, and 
its well known manifestation is current induced mag- 
netic switching (CIMS) 3 -^ as well as generation of mi- 
crowave oscillations 5 -* 6 -. Properties of STT are related 
to the sample design and material parameters. This 
also implies that CIMS and current-pcrpendicular-to- 
plane giant magnetoresistance (CPP-GMR) phenomena 
are correlated^ and depend on the same structural pa- 
rameters. By considering two well-defined conduction 
spin-channels, Valet and Fert£ incorporated most of these 
parameters into CPP-GMR model based on spin diffu- 
sion transport equations. Generalization of Valet-Fert 
approach^ includes STT and provides an unified descrip- 
tion of STT and CPP-GMR in the diffusive transport 
limit. For symmetric spin valves, with fixed and free 
layer made of the same material, such as Co/Cu/Co, 
this model predicts standard behavior of STT, which 
does not vary qualitatively from Slonczewski's result ob- 
tained in the ballistic transport limiti In this case, 
current of one orientation drives switching to antipar- 
allel configuration while opposite current stabilizes par- 
allel stated. In presence of applied fields higher than 
the coercive field, the generation of microwave oscilla- 
tions is possible^. Similar behavior has been recently 
observed also in Py/Cu/Py (Py = Permalloy), both 
exp eriment all y 1 1 1 1 2 and theoretically^ 3 -. Moreover, stan- 
dard STT associated precession was also reported at 



high perpendicular magnetic field in an asymmetric spin 
valvei^, where the magnetic layers are made of different 
materials (Co(40 nm)/Cu(10 nm)/Py(3 nm)). 

However in such asymmetric structures, provided 
that spin asymmetries and spin diffusion lengths dif- 
fer markedly and the thickness obeys certain conditions, 
qualitatively different situation may arise. The STT van- 
ishes and changes sign in a certain noncollinear magnetic 
configuration (wavy-like STT) due to the appearance of 
an inverse spin accumulation in the nonmagnetic spacer. 
Therefore, current flowing in one direction destabilizes 
both collinear magnetic configurations, whereas the op- 
posite current stabilizes both of them. The first case is of 
particular interest as it leads to excitation of stationary 
oscillation modes in absence of external magnetic field. 

In this paper we present a systematic study of dynamic 
response of a magnetic film to such wavy-like STT. The 
asymmetric pillar under study is a spin valve consisting 
of a fixed layer, nonmagnetic spacer and a free magnetic 
layer, Co(8 nm)/Cu(10 nm)/Py(8 nm) respectively, with 
the elliptical cross-section of 155 x 100 nm 2 . The polar- 
ization of the fixed layer is assumed to be along the ellipse 
major axis. Recently, such asymmetric structures have 
been investigated theoretically ^ 15 ' 16 as well as experi- 
mentally at low and zero applied field a 17 : 18 . To authors' 
best knowledge, no micromagnetic analysis of such struc- 
tures in the diffusive transport limit has been performed 
so far. Moreover, the reported macrospin (referred to 
in the following also as single domain) study does not 
describe correctly the dynamics at low applied magnetic 
fields^. Here, starting with the single domain approxi- 
mation, and extending study to the full micromagnetic 
model, we explain the origin of out-of-plane precession 
(OPP). In particular, we show that only full micromag- 
netic model can successfully reproduce magnetization dy- 
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namics at low applied magnetic field. 

The paper is organized as follows. Section [TT1 describes 
briefly the torque calculations in the diffusive transport 
limit. The methodology of simulations is presented in 
section IIII1 The results of numerical study and their 
discussion are to be found in section IIV( whereas final 
conclusions in section Ivl 



II. SPIN TORQUE IN AN ASYMMETRIC 
PILLAR 

As mentioned in the introduction, structure of a pillar 
determines the dependence of STT on the angle between 
magnetization vectors. Generally, STT consists of two 
components, r = tq + t v , where Tg is the in-plane (IP) 
component, while t v is the out-of-plane (OP) one. These 
two components can be written as^ 

tq = —aj m x (m x p) , (la) 
= bj m x p , (lb) 

where m denotes the normalized (unit) vector along the 
magnetization of the free layer, p is the normalized mag- 
netization of the fixed layer, and j is the current density. 
The prefactors a and b are independent of current j, but 
they generally depend on the angle 9 between m and p. 
These parameters have been computed from the diffusive 
transport model^. First, from the boundary conditions 
for spin current and spin accumulation^ we calculate the 
spin current. The torque is then calculated from the nor- 
mal component of the spin current in the nonmagnetic 
film at the interface with the magnetic layer. Most of the 
parameters used in this description, like interface and 
bulk spin asymmetry coefficients, interface resistances, 
layer resistivities^, and spin diffusion lengths^ are pro- 
vided by the corresponding CPP-GMR experiments. The 
two remaining parameters, i.e. the real and imaginary 
part of the mixing conductance^ have been extracted 
from spin current interface transmission calculations^. 
The variation of normalized (to Hj/\e\) STT with the an- 
gle 9 for symmetric and asymmetric spin valve is shown 
in Fig. m 

In the system under considerations (Fig. [T|), the OP 
torque is roughly 2 orders of magnitude smaller than 
the IP component, which means that the latter will 
not markedly influence the magnetization dynamics. 
Comparing the angular dependence for the standard 
(Fig. QJ,) pillar structure with that for the nonstandard 
one (Fig. [TJd), the uniqueness of the latter is clearly vis- 
ible - at some critical angle 9 the IP component of the 
torque vanishes. This gives rise to interesting dynamics 
at zero and low magnetic fields. Above certain threshold 
current both collinear states of the magnetization are un- 
stable for one current orientation, and the only solution 
of the Landau-Lifshitz-Gilbert (LLG) equation is then 
the steady state precession or a noncollinear static mag- 
netization stated. 
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FIG. 1: Angular dependence of the spin trans- 
fer torque normalized to hj/\e\ for (a) symmetric 
Py(8nm)/Cu(10nm)/Py(8nm) and (b) asymmetric 
Co(8nm) /Cu(10nm)/Py(8nm) pillar. 

III. SIMULATION METHODOLOGY 

Magnetic dynamics of the free layer is described by the 
Landau-Lifshitz-Gilbert (LLG) equation with the STT 
included. This equation, when written as the equation 
for time evolution of the unit vector m, takes the form 



1 + a 2 dm 

-— — = -m x h cf f - am x (m x h off ) 

70 M s dt 

+ ^ M i d I am x (m x p) - 6m x p] (2) 

where a = a — ab, b = b + aa, h e ff is the effective field 
normalized to the saturation magnetization Ms of the 
free layer, a is the damping constant, d is the thickness 
of the free layer, and 70 = |'y|A t o with 7 being the gyro- 
magnetic ratio and fiQ standing for the magnetic vacuum 
permeability. Since a is a small parameter, in numerical 
simulations we assumed a — a and b = b. 

In our convention the current flowing from the fixed 
towards the thin (free) layer defines the positive cur- 
rent direction, j > 0. Note, that in this convention 
the positive current triggers switching to the antiparal- 
lel configuration, whereas negative current stabilizes par- 
allel alignment in standard spin valves. On the other 
hand, negative current in nonstandard spin valves sup- 
ports oscillatory regime, while positive current stabilizes 
both collinear configurations. 

Only the dynamics of the free layer is resolved, while 
magnetization of the fixed layer is assumed to remain uni- 
form in the film plane. Apart for this, we have neglected 
all thermal effects. In the macrospin analysis, the ef- 
fective field is assumed to include the self-magnctostatic 
terrn^, the uniaxial anisotropy field, and the external 
magnetic field. The uniform magnetostatic coupling field 
(with the fixed layer) can also be taken into account. In 
the micromagnetic study, on the other hand, the mag- 
netic field is modified with respect to macrospin case. 
First, the self-magnctostatic term is computed using 
the Fast Fourier Transform technique assuming that the 
magnetization is uniform in each computational cell2£. 
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Second, six-neighbor dot product representation is used 
to compute the exchange field. Unless stated differently 
in the text, the following values of the relevant param- 
eters have been chosen for the systematic study: the 
anisotropy constant K u = 3.46 x 10 3 J/m 3 , exchange con- 
stant A = 1.3 x 10~ n J/m, damping constant a = 0.01, 
M s = 6.9 x 10 5 A/m, and M sfixod = 1.4 x 10 6 A/m. As 
follows from this choice for permalloy, the corresponding 
exchange length Z cxc h = y / 2A/ ^ M'g equals 6.6 nm, and 
therefore we choose 5x5x4 nm 3 discretization mesh. 
Further refinement of the cell size in the z-direction does 
not lead to any substantial difference in system frequency 
response. Fourth-order Runge-Kutta scheme was em- 
ployed for the time integration of Eq. [21 and the stability 
analysis was carried out for the chosen mesh to assign 
appropriate time integration step. 



IV. RESULTS AND DISCUSSION 



glected. Since the anisotropy of Py is also small, one 
concludes that it is mainly the self-magnetostatic field 
that drives the system dynamics. In the following we 
consider the role of initial magnetic state of the system, 
and begin with the parallel magnetic configuration. 
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FIG. 2: Geometry of a pillar with (a) an extended and (b) 
etched fixed layer. The interlayer coupling field in (a) can be 
neglected. 



Before comparing results obtained from macrospin ap- 
proximation with those from full micromagnetic study, 
one should consider that the concept of single domain 
magnetic particle in many situations is not justified. We 
shall explain how the difference between the results of 
macrospin and micromagnetic analysis arises from the 
appearance of an inhomogencous magnetization and a fi- 
nite exchange field. The exchange energy density of a 
closed M(r) configuration increases as the particle size 
decreases, and this could, in principle, justify macrospin 
approach below certain critical size of the system, even 
though spin torque makes the estimation of this criti- 
cal size difficult. However, as reported in [28[, steady 
state precession of a thin square nanoelement exhibits 
complicated transition from quasi-macrospin to chaotic 
behavior already at the size of 30 nm, which invalidates 
single domain approximation for most of experimentally 
studied systems. Moreover, on the basis of micromag- 
netic analysis Berkov and Gorn^S identified some arti- 
facts of macrospin model in the ballistic transport limit. 
These artifacts might cause misleading interpretation of 
the origin of some observed phenomenal. Aware of that, 
we cannot however deny that macrospin dynamics is very 
rich and in many cases can serve as a tool for the basic 
understanding of physics behind. 

Since the precessional states appear for negative cur- 
rent (according to the definition introduced in section 
III) , we limit the following considerations to negative cur- 
rent. For sake of simplicity, the current / denotes the 
absolute value of the negative current. 



A. Extended geometry 

We start our numerical study with a pillar struc- 
ture having an extended fixed layer, like the one shown 
schematically in Fig. [2k- In this case, the influence of 
interlayer magnetostatic coupling field (ICF) can be ne- 



1. Initial P state 

Assume the system is initially in the parallel configu- 
ration and increase the current step-like (with the step 
AI) in order to identify possible oscillatory regimes. As 
the initial state for simulation at a given current / we 
assume the state the systems arrived at in the preced- 
ing step of simulations, i.e. for current / — AI. Sim- 
ulation results from both macrospin approximation and 
micromagnetic model reveal that an in-plane precession 
(IPP) is supported at zero applied field, which is con- 
sistent with previous study flMlril ]. The associated fre- 
quency redshift with current is a well known phenomenon 
and is attributed to the increase of the oscillation am- 
plitude. Interestingly the redshift in the micromagnetic 
model turns out to be much slower than in the macrospin 
approximation, which is opposite to the results obtained 
in the ballistic transport limit^l. We shall explain this 
effect later on. 

As follows from Fig. [31 the precessional states in the 
macrospin model disappear above a certain 'cut-off' cur- 
rent (1=35 mA), where a stable static 'spin-up' (normal 
to the film plane) magnetization state is formed. Ex- 
istence of this static state has been recently reported in 
mum, where it was found that the LLG equation has 
just two possible solutions - self-sustained precession or 
a stable static state (above the 'cut-off' current). The 
appearance of the latter can be explained with the help 
of Fig. [4] As the current is applied and the system is ini- 
tially in the P state, the STT counterbalances the damp- 
ing, and steady state IPP+ is obtained in region I. When 
the current is increased, the oscillation amplitude in- 
creases and therefore the critical angle (marked as cross), 
at which the torque vanishes, is reached. The magneti- 
zation starts then aligning along the effective magnetic 
field, and finally the stable 'spin-up' state is reached. 
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No such static state has been reported in micromagnctic 
simulations. Assume now that the system in macrospin 
simulations is in the static 'spin-up' state and then the 
current is decreased (region II). The self-magnetostatic 
field in the static point together with the STT may trig- 
ger the OPP in some current range, as shown in Ref. 
[16l | . Again, no OPP appears in the micromagnetic sim- 
ulations. In general, when the angle between magnetiza- 
tions of the free and fixed layer is lower than the critical 
angle (region I), only IPP is observed. As the critical an- 
gle is crossed (region II) the OPP can be triggered in the 



I5II32I ]. However, we 



macrospin model, as shown in Refs. 
point again that the 'spin-up' state is not reached in mi- 
cromagnetic simulations because of the inhomogeneous 
character of the magnetization. Therefore, even though 
single domain model predicts OPP, which was also re- 
ported experimentally^, the interpretation of its origin 
requires further considerations. 
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FIG. 4: The torque normalized to frj/\e\ for negative current 
density (j < 0, compare Fig.l) acting on the free layer. As 
the current is increased, the frequency redshift due to the am- 
plitude increase is observed (angle increase, region I). When 
the critical angle is reached (cross between the regions I and 
II), the torque vanishes and the magnetization tends to align 
along the effective field driving system out of the critical angle 
towards the stable static 'spin-up' state. 
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FIG. 3: Frequency vs. current behavior in absence of external 
field in the macrospin model - full triangles, and in micromag- 
netic model - open squares. This convention is kept through- 
out this paper: open symbols correspond to results of micro- 
magnetic study and full symbols to those in the macrospin 
model. A 'cut-off' current is observed at 1=35 mA in the sin- 
gle domain model. The simulations have been performed for 
the step- wise increasing current, and for parallel initial con- 
figuration of the system. Initial state for a current I + AI 
was taken as the state at which the system arrived at in the 
simulations for the current /. 



2. Initial AP state 

As we have seen above, two factors play crucial role in 
the excitation of precessional states in the system. These 



are the magnetization state which imposes the initial self- 
magnetostatic field, and the initial angle between the 
magnetic layers, which determines the initial torque sign 
and strength. The IPP is supported in region I (Fig. 
whereas in region II the static state has been found in 
the macrospin approximation so far. 

In order to investigate modes supported in region III, 
we assume AP configuration as the initial state. As the 
current is increased, micromagnetic simulations give a 
fast red-shifting branch corresponding to IPP_ (preces- 
sion around the —x axis) - marked as 1 in Fig. O (open 
squares). Fast amplitude increase with increasing current 
soon leads to switching towards the P state and damped 
oscillations stabilizing P state are observed (range 2). 
However, as certain threshold current is reached (P is 
now the initial state) , the STT (region I) counterbalances 
the damping and the second red-shifting branch, IPP+, 
is observed (range 3). Micromagnctically, the only differ- 
ence between starting with initial P or AP state (open 
circles and open squares respectively) is the appearance 
of first fast red-shifting branch. 

Analogous comparison (see full circles for initial P state 
and full squares for initial AP states, Fig. [5]) in the 
macrospin approximation leads to similar conclusions, 
but now an additional mode is visible. When current is 
high enough and assuming the AP state as the simulation 
is initialized, the system might be forced to move into re- 
gion II leading to the appearance of OPP marked as 4 
(full squares). Earlier in this section it was shown that 
after crossing from region I through the critical angle into 
region II, the static state was observed. Here, however, 
dynamics in region II is forced by the initial configura- 
tion and therefore the OPP can be observed. In other 
words, in the single domain approximation the region I 
supports IPP+, region II supports static state or OPP 
(depending on the proceeding configuration), whereas in 
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region III the IPP_ mode can be observed. This result 
is consistent with the one reported in [l5l |. 



12 
10 



o 

o 
a 

V 

3 
a* 

o 4- n 



□ m m = (-1,0,0) 
O m'° = ( 1,0,0) 




10 



20 30 
current[mA] 



40 



FIG. 5: The influence of the initial magnetization state on the 
dynamic behavior of the system and the corresponding spots 
in torque diagram (inset). Circles refer to P, whereas squares 
to AP initial state. At low currents, IPP_ mode (close to 
the AP direction) is found (marked as 1). As the current 
increases, the angle of the orbit increases as well, and (m x ) 
switches towards the P state. No self-sustained oscillations 
are observed in the area 2. Further current increase leads to 
the onset of the second red-shifting branch, IPP+, marked 
with 3. Additionally at a certain threshold current, the com- 
bined effect of self-magnetostatic field of the AP initial state 
and negative torque sign can trigger OPP in the macrospin 
model, marked as 4. 



B. Etched geometry 

In the extended structures (Fig. [2^) as discussed above, 
the interlayer coupling field (ICF) could be neglected, 
and the initial self-magnetostatic field together with the 
STT determined the induced dynamics (region I, II 
or III). However in the etched geometry, schematically 
shown in Fig. ^Bp, the ICF can no longer be neglected. 
We have calculated this field micromagnetically, and by 
neglecting large OP edge values we have estimated -36mT 
IP as an average ICF. This is a significant contribution 
and therefore dynamics different from that obtained for 
extended structures is expected. 



1. Initial P state 

The dynamical response reveals now some new inter- 
esting features. The macrospin approximation with the 
initial P state leads to the steady OPP, which appears at 
10.4 mA when the current is increased, as indicated in 
Fig. [5] (corresponding range for observed mode marked 
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FIG. 6: The influence of current increase and decrease in the 
etched pillar geometry on the dynamical response of the sys- 
tem with P initial state. The corresponding torque diagram 
is shown in the inset. Note asymmetry in the macrospin fre- 
quency response - OPP with current increase and transition 
from OPP to IPP with current decrease are observed. This 
hysteretic behavior originates from the asymmetric torque an- 
gular variation and the existence of the critical angle making 
dynamic transition between region I and II prohibited. 



with dotted lines). As the threshold current for OPP 
is reached, the blue-shifting branch appears in range 2. 
This OPP is not preceded by any IPP oscillations - range 
1, because the additional contribution from ICF places 
the system directly in region II. As the frequency of OPP 
oscillations increases with increasing current, the corre- 
sponding amplitude decreases and the angle between the 
magnetic moments of both layers approaches the critical 
angle. When this angle is reached, the static state dis- 
cussed in the preceding subsection is observed (point 3). 
Now we start to decrease current. The system is initially 
in the static point and dynamic range is marked with 
solid lines in Fig. [6l The OPP appears then in the range 
1, and the amplitude and 6 increase as the current de- 
creases. The torque minimum is then passed and the sys- 
tem moves to the region III, which results in IPP_ mode 
in range 2, where the amplitude decrease (with decreas- 
ing current) results effectively in the frequency rcdshift 
with current. The asymmetry in the macrospin frequency 
response to increasing and decreasing current is clearly 
a consequence of the torque asymmetry, the existence of 
a critical angle, and irreversibility of the transition from 
region I to region II. 

Dynamics in the micromagnetic model is simpler. Due 
to the effect of ICF, the system directly switches to 
the region III, and only one red-shifting branch IPP_ 
is observed. This mode is qualitatively equivalent to the 
macrospin one in range 2 observed for decreasing current. 
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FIG. 7: The same as in Fig. [6] but for the AP initial state. 
Now no frequency hysteretic behavior with current is observed 
in the macrospin model, and the dynamics evolves reversibly 
between the regions II and III (inset). Micromagnetically only 
the IPP mode in the region III is observed. 



netically due to the inhomogcnous nature of the mag- 
netization (finite exchange field) in the model. There- 
fore, one should expect that the appearance of OPP in 
micromagentic model is hindered by the underestimated 
exchange field, and that increase of the exchange con- 
stant should lead to a convergence of both models. Using 
the bulk exchange constant for thin films might cause a 
great underestimation of exchange fields in these struc- 
tures. Larger values of exchange constants, as compared 
to the standard bulk ones, have been reported in Py 
dots^ and thin films24. Therefore, we have investigated 
magnetization dynamics in the extended geometry for 
the following values of the exchange constant: 0.75A , 
2.5A , 3A , 4Aq, 10 Aq. The frequency-current behavior 
with A = 3Aq is compared to the results of macrospin 
model in Fig. [51 Clearly increasing the exchange constant 
changes the slope of micromagnetic frequency rcdshift 
towards macrospin results. We conclude that the finite 
exchange energy favoring inhomogenous magnetization 
state causes this slope difference retarding the dynamics. 



2. Initial AP state 

As discussed above, starting from the initial P state in- 
troduced in the macrospin approximation hysteretic de- 
pendence in frequency vs. current response. Clearly it is 
due to the asymmetric shape of the torque angular de- 
pendence. This effect, however, disappears when the ini- 
tial state is AP. When the current is increased (dotted in 
Fig. [7]), the system is directly placed in region III support- 
ing red-shifting branch (IPP_) in range 1 (inset Fig. EJ). 
As the amplitude increases with increasing current, tran- 
sition to range 2, where OPP oscillations are triggered, 
is observed. In the torque diagram this is equivalent to 
the transition from region III over the torque minimum 
to region II. Since the critical angle is not crossed, this 
transition stays reversible and no hysteretic behavior in 
the frequency response is observed. 

As before, micromagnetically no OPP was found. This 
indicates that the system supports stable oscillations only 
in region III. 
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FIG. 8: The dynamic response of the system for two different 
values of the exchange constant, compared to the macrospin 
results. Micromagnetically faster redshift is observed for 
larger exchange constant as it favors more homogenous mag- 
netization and accelerates the dynamics. 



C. Influence of the exchange field 

1. The extended geometry 

An open question is why the OPP modes obtained in 
the macrospin model and also reported experimentally 18 , 
have not been found in micromagnetic simulations. As 
has been shown, to observe the OPP-associated blucshift , 
one has to force the system dynamics in region II. More- 
over, we have learnt that the appearance of the macrospin 
'spin-up' state was a consequence of system crossing over 
the critical angle. This has not been reached micromag- 



2. Etched geometry 

So far in this geometry (ICF included) we have found 
that as the current was increased micromagnetically in- 
duced dynamics in region III imposed frequency rcdshift 
and the macrospin dynamics in region II supported OPP- 
associated blueshift. Still we want to check whether mi- 
cromagnetic dynamics can be shifted to region II by in- 
creasing the exchange constant. The micromagnetic tem- 
poral evolution of the averaged magnetization component 
at 1=13 mA for A — Aq and A = 2.5Ao results in dif- 
ferent orbits, Fig. [9fi and Fig. O), respectively. Clearly 
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FIG. 9: Visualization of the exchange constant influence on 
the magnetization dynamics in the etched geometry at 1=13 
mA. Magnetization orbit for (a) A = Ao typical for region III, 
and (b) A = 2.5Ao approaching the border between regions 
III and II 
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in the first case the ICF places the system in region III 
forcing IPP_ dynamics. However, as the exchange field is 
increased, which favors uniform magnetization, an open 
clamshell orbit is formed (Fig. [QJd) shifting the dynamics 
towards the border between region III and II. One should 
note that the cross-over between the regions is impossible 
in this geometry since the ICF has a significant contri- 
bution to the effective field and hinders the appearance 
of OPP. 



FIG. 10: Comparison of frequency vs. current behavior at 
OmT between macrospin model and micromagnetic model 
with A = 3Ao. As described in the text, in the micromag- 
netic model the torque strength has been scaled by a factor 
of 0.5 in order to counterbalance the effect of its inhomoge- 
neous nature. The threshold for OPP predicted here by both 
models coincides with the experimental one [is| . The mag- 
nitude of the frequency jump associated with the transition 
IPP-OPP in the macrospin approximation does not fit the ex- 
perimental values and micromagnetic approach proves to be 
more accurate. 



D. Comparison to experimental data 



Still no OPP (supported in the dynamic region II) has 
been predicted micromagnctically, although such modes 
were reported experimentally^ in extended structure. 
However previous paragraphs have given some important 
clues. We have learnt about the importance of the under- 
estimation of the exchange field. Therefore, we assume 
A = 3Aq for further study. Secondly, as the magneti- 
zation always stays inhomogeneous to some extend, the 
toruqe calculated locally (cell by cell) inherits this inho- 
mogencity and we shall scale the torque strength by a 
factor of 0.5 to counterbalance this effect. Thirdly, mi- 
cromagnctically the transition from dynamic region I to 
II was impossible. Thus in order to observe OPP one has 
to force dynamics in region II by forcing the transition 
from region III to II (i.e. imposing AP initial state). As 
presented in Fig. 1 101 indeed under all above mentioned 
assumptions both models converge. Micromagnetic dy- 
namics is forced first in region III supporting IPP_, then 
switching towards P state takes place (range of current 
where no sustained oscillations are observed, as discussed 
before), and then IPP + branch (region I) is triggered. 
At a certain threshold, however, the dynamics in region 
II supporting blucshift can be obtained. The threshold 
of this OPP coincides with the experimental one from 
[l8j]. The magnitude of the frequency jump associated 



with the transition IPP-OPP in the macrospin approxi- 
mation does not fit to the experimental values, and mi- 
cromagnetic approach proves to be more accurate. The 
fact that experimentally OPP was reported by starting 
from P state (opposite to our results) means that tran- 
sition between region I and II prohibited micromagnct- 
ically is experimentally possible due to thermal activa- 
tion. As our simulations neglect the effect of thermal 
fluctuations, the dynamic region II can be only reached 
by transition from region III. One should note that not 
all regimes (dynamics in all regions I, II, III) predicted 
by the simulations for the extended geometry were ob- 
served experimentally. Low angle IPP does not provide 
enough output power to be measured via GMR effect. 
Therefore in order to conduct meaningful comparison we 
have concentrated on the OPP regime (region II) , which 
was both predicted numerically and observed experimen- 
tally. In the absence of external field satisfactory quali- 
tative agreement has been reached (Fig. [Til compare to 
Fig. 6a in Ref. [HI]). Not only the threshold current 
(Ith,sim = 11 mA compared to Ith.exp = 10 mA) but also 
the agility (0.6 GHz/mA and 0.7 GHz/mA respectively) 
are in good agreement. The remaining quantitative dif- 
ference in frequency values is a consequence of uncertain 
estimation of the factors entering micromagnetic model, 
like saturation magnetization and/or damping. More- 
over, the dynamics at -5mT IP field reveals that OPP 
threshold current is smaller with respect to OmT case, 
which is again consistent with the experimental results. 
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of low output power) so obviously the experimental cut 
off current refers to the threshold current for IPP reap- 
pearance in the model. Note that micromagnetically the 
main mode (supported over largest range of currents) in 
case of OmT was the IPP+, whereas even low applied 
field of -5mT IP forced the dynamics in region III and 
therefore the IPP_ was observed as the main mode. In 
other words in the absence of external field the increased 
exchange constant enabled direct dynamics in region II 
(and associated linear frequency vs. current slope) and 
additional field forced the transition from region III to 
region II resulting in the appearance of the saturation 
regime. Clearly the OPP in both cases is preceded by 
different dynamics. 



FIG. 11: Results of micromagnetic model at zero and low 
applied field. The dynamic response at OmT and -5mT ap- 
plied in-plane. In absence of external field IPP+ transition 
into OPP at 10.5mA and back to IPP+ above 11.5mA is ob- 
served. At -5mT the transitions from IPP_ to OPP and back 
to IPP_ are observed at 9.5mA and 12mA respectively. Both, 
blueshift linear at OmT and nonlinear at -5mT predicted mi- 
cromagnetically, are consistent with the experimental results 
[l|| . Macrospin ceases to show the saturation regime observed 
experimentally at -5mT but full micromagnetic study cor- 
rectly predict this feature. 

Interestingly, both approximately linear at OmT and non- 
linear at -5mT behavior of the frequency as a function of 
current, are well reproduced in frames of micromagnetic 
model. Note that this feature was not reported in single 
domain model. Furthermore, micromagnetic model pre- 
dicts experimentally observed saturation at -5mT, i.e. in 
a certain current range the frequency stays relatively con- 
stant. In the model it is associated with the large angle 
orbit stabilization around the torque minimum, i.e. the 
system approaches the border of dynamic regions II and 
III, and the torque shape becomes flat around its mini- 
mum. Furthermore, current increase leads to transition 
between OPP (region II) and IPP_ (region III) and reap- 
pearance of the clamshell orbit. However, since neither 
was the IPP reported in the experiment prior to the ap- 
pearance of OPP, as predicted micromagnetically, nor it 
could have been detected following OPP regime (because 



V. CONCLUSIONS 

The results emerging from the macrospin and micro- 
magnetic model very often do not converge. Lack of this 
convergence is attributed to a couple of factors. It was 
shown that the discrepancies arise from the inhomoge- 
nous magnetization state and finite exchange field. We 
concluded that the OPP in the macrospin model might 
appear as a consequence of the 'spin-up' static state state, 
which seems to be characteristic of the macrospin model. 
On the other hand, the absence of OPP in the micromag- 
netic model was identified as a combined consequence of 
the underestimation of the exchange constant and the 
role of the initial self-magnetostatic field. By setting the 
initial state and the exchange constant favoring the ap- 
pearance of the OPP, a good qualitative agreement was 
reached between the predictions of both models. Thus, 
only micromagnetic model has proven to predict correctly 
dynamics reported experimentally. 
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